Ground-based detections of thermal emission from CoRoT-lb and 

WASP-12b 



Ming Zhao 1 ' 2 , John D. Monnier 3 , Mark R. Swain 1 , Travis Barman 4 , Sasha Hinkley 5 ' 1 

ming . zhao@ j pi . nasa . gov 

ABSTRACT 

We report a new detection of the if-band thermal emission of CoRoT-lb and 
two confirmation detections of the Ks-band thermal emission of WASP- 12b at 
secondary eclipses. The ff-band measurement of CoRoT-lb shows an eclipse 
depth of 0.145%±0.049% with a 3-a percentile between 0.033% - 0.235%. This 
depth is consistent with the previous conclusions that the planet has an isother- 
mal region with inefficient heat transport from dayside to nightside, and has a 
dayside thermal inversion layer at high altitude. The two Ks band detections of 
WASP-12b show a joint eclipse depth of 0.299%±0.065%. This result agrees with 
the measurement of Croll & collaborators, providing independent confirmation 
of their measurement. The repeatability of the WASP-12b measurements also 
validates our data analysis method. Our measurements, in addition to a number 
of previous results made with other telescopes, demonstrate that ground-based 
observations are becoming widely available for characterization of atmospheres 
of hot Jupiters. 
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Introduction 



Detecting thermal emission from transiting planets at secondary eclipses is a power- 
ful technique to study their atmospheres without spatially resolving them from their host 
stars. Recently, ground-based observations have emerged as a powerful tool, in addition 
to sp ace-based observations, to characterize the a t mospheres of transiting hot Jupiters 
(e.gJsing fc Lopez-Morale] l2009l : IRogers et allbood ; lLopez-Morales et alJlioiot ICroll et al. 
2010al etc.). These ground-based photometry observations mostly provide measurements at 
the z' , J,H,& K band, with the near-IR (J, H, K) measurements probing dee per and higher 
press ure layers of planetary atmospheres than the Spitzer measurements (IBurrows et al. 
20081). Since the b ulk of energy from hot Jupiters emerges from the near-IR between 1-3 
fim (IBarmanl 120081 ). ground-based observations have the potential to provide an important 
probe of these atmospheres. 

CoRoT-lb and WASP- 12b are two Very Hot Jupiters that are well suited for ground- 
based secondary eclipse detections. CoRo T-lb 



has emer gent flux detected in both optical ([Snellen et al. 



Barge et al. 2008h is the first pla net that 



20091 ; lAlonso et all 120091 ) and the 



near-IR (IGillon et al.l 120091 ; IRogers et al.l 120091 ) . Recent ly, its thermal emis sion at 3.5 and 
4.6 \im has also been measured by the Spitzer telescope (IDeming et al.1120111 ). These studies 
provided multi-wavelength constraints of CoRoT-lb's atmosphere, suggesting it has a tem- 
perature of ~ 2310K and possesses a therma l inversion layer at h igh altitude, consistent with 



the pM-class classification of hot Jupiters (iFortnev et al 



shown to have an inflated radius of !A5Rj up (IGillon et al. 



2008). CoRoT-lb has also been 



20091 ). which could be caused by 



tidal heating. However, its orbital ecce ntricity has been shown to be consistent with zero 
(IRogers et aL||2009j ; IDeming et al.1120111 ) . thus not providing any support for the theory that 
CoRoT-lb's inflated radius is due to the effects of tidal heating. 

Despite the extensive studies of CoRo T-lb, current atmospheric models stil l under- 
predict its Ks and 2.09/im measurements (IGillon et al.l 1200 9; Ro gers et al.l 120091) . and a 



blackbody model provides the best fit to the measurements (IDeming et al 



20111 ). Additional 



measurements at other bands such as H are desirable to further constrain CoRoT-lb's models 
and help us better understand its atmospheric characteristics. 



WASP- 12b is one of the hottest and most inflated planets discovered to date (IHebb et al. 



2009 ). Powerful irradiation from its host star heats the atmosphere temperature above 2500K 
(IMadhusudhan et al.l 120111 ). resulting in deep eclipse depths that are favorable for ground- 
based secondary eclipse detections. Currently, thermal emission of W ASP- 12b has been de 



tected in the z', J , H, K bands and the Spitz er 3.5/im and 4.6/im bands (lLopez-Morales et al. 



2010 



Croll et al 



2011 



Campo et al.l 1201 ll ). Based on these observations, recent models of 



Madhusudhan et al.l (120111 ) suggest that the atmosphere of WASP- 12b is extremely carbon- 
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rich (specifically the C/0 ratio of WASP-f2b is >1 at 3-a significance). However, it lacks a 
prominent thermal inversion layer at photospheri c depths predicted fo r very-hot Jupiters and 
has very efficient day-night heat redistribution (IFortney et al.l 120081 ). These characteristics 
motivate new planetary inte rior models and present ch allenges to theoretical classifications 
of hot- Jupiter atmospheres (IMadhusudhan et al.l 1201 fl ). 



Here we report a new detection of CoRoT-lb's thermal emission at the H band with 
the Palomar 200in telescope, and two confirmation detections of WASP-12b's Ks band 
emission with the Michigan-Dartmouth-MIT (MDM) 2.4m telescope. In §5] we present our 
observations and data reduction procedure. We discuss our data analysis process and the 
results of CoRoT-lb and WASP-12b in §3j In §4]we compare the eclipse depth of CoRoT-lb 
with existing models, and finally we summarize our results in §5J 



2. Observations and data reduction 
2.1. CoRoT-1 

The observation of CoRoT-lb was conducted in the H band with the WIRC instrument 



flWilson et al.l 120031 ) on Palomar 200-in Hale telescope on UT 2011 January 29. The WIRC 
camera has a 2048 x 2048 Hawaii-II HgCdTe detector with a pixel scale of 0.2487"/pixel 
and a wide field of view (FOV) of 8.7' x 8.7'. The observation started about 10.5 minutes 
before the ingress and ended about 102.5 minutes after the egress, lasting for about 255 
minutes in total. Each image was taken with 30sec exposure and 1-fowler sampling. A total 
of 376 images were obtained. The duty cycle of the observation is 73%. The average seeing 
of the night was ~ 0.8". We stayed on the target for the whole period without dithering 
to minimize instrument systematics. The telescope was slightly defocused to a FWHM of 
~ 1.5" - 2" to keep the counts well within the linearity regime and to mitigate any potential 
intra-pixel variations. We adjusted the telescope occasionally throughout the observation to 
maintain the defocus. Because WIRC has no dedicated guider, the telescope pointing drift 
was manually corrected during the observation. 

For the data reduction process, we first subtracted all images with corresponding av- 
eraged dark frames. Sky flats were normalized and averaged to get a master flat field. We 
create a bad pixel mask with the master flat and dark frames. The bad pixels in each im- 
age were interpolated with cubic splines based on adjacent flat-fielded pixels. After these 
steps, stars within the flux range of 0.25 to 1.5 times of that of CoRoT-1 are selected as 
references. Stars with higher fluxes are beyond the linearity regime of the detector and thus 
are excluded. Stars with lower fluxes have insufficient signal-to-noise and are excluded as 
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well. In addition, 4 stars within this flux range are excluded due to their excessive flux 
fluctuations compared to other stars. The selection leads to 31 well separated and evenly 
distributed reference stars in the FOV for flux calibration. Due to the correlation of stellar 
flux variations with their centroid positions on the detector (see §3.ip . more reference stars 
are preferred to fewer stars in order to average out their correlation with centroid positions 
on the detector. We calculated the centroids of all stars in each image using a center-of-mass 
calculation, since it provided the smallest scatters of their relative positions. The time series 
of CoRoT-l's centroid was determined by averaging the relative positions of all reference 
stars after correcting for their relative distances. The resulting 1-a precision of the centroid 
determination is ~ 0.3 pixels. 

Aperture photometry was performed on CoRoT-1 and the reference stars following the 
IDL routines of DAOPHOT. The extracted fluxes of each star are normalized to the median 
of the time series. The median^ of the 31 reference time series is then taken as the final 
reference light curve, which is used later to normalize the flux of CoRoT-1 to correct for 
the common-mode systematics such as variations of atmospheric transmission, change of 
seeing and airmass, etc. We applied 48 different aperture sizes with a step of 0.5 pixels, 
and determined that an aperture with a radius of 9.5 pixels (19-pixel diameter) gives the 
smallest out-of-eclipse and in-eclipse scatters for the normalized CoRoT-1 data; this is taken 
as the final photometry aperture for all stars in every image. Apertures within radii of 
9.5 ± 1.5 pixels show consistent eclipse depths in later analysis, while apertures with larger 
than 1.5-pixel differences start to show excessive systematic noises. A sky annulus with 35- 
pixel inner radius and 30-pixel width was used for background estimation. The median value 
of the sky annulus was then used as the final sky background for subtraction. We have also 
explored different annulus ranges and sizes, and found consistent results. The top two panels 
of Figure [1] show the reduced fluxes of all 32 stars and the final normalized flux of CoRoT-1. 
The relative centroid changes of CoRoT-1 are shown in th e two bottom panels. The UTC 



timestamp of each image was converted to BJDutc fir st (Eastman et al.ll2010l ). and then 



converted to orbital phases based on the ephemerides of ( Bean 20091 ). i.e., period=l. 5089656 
days, and transit epoch T (BJD UTC ) =2454159.452879. 



^^Due to the presence of outliers in the light curves of some reference stars, as can be seen in Figure 1, we 
use median here as a more robust estimator. 
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Fig. 1. — Reduced fluxes and relative positions of CoRoT-l's centroid. The first panel 
shows the normalized fluxes of all 32 stars, including CoRoT-1 (black dotted line in the top 
panel). The second panel shows the final CoRoT-1 data normalized with the median of all 
reference fluxes. The two bottom panels show the relative positions of CoRoT-l's centroid. 
The centroid varies by >10 pixels during the observation. The blue line in the second panel 
shows the best-fit de-correlation model combined with a best-fit light curve (de-correlated 
with X positions only, see §3. lj) . The flux variations in the top panel of Figure [D are due 
to systematic effects such as airmass, seeing variations, and fluctuations of the atmospheric 
transmission, etc. 



2.2. WASP-12 



The observations of WASP-12b were conducted on two nights (UT 2010 November 26 
& 27) at the Ks band with the TIFKAM imager on the MDM 2.4m Hiltner telescope. 
TIFKAM has a 1024x1024 HAWAII-1R HgCdTe detector. We used the f/7.5 imaging stop 
with a pixel scale of 0.2"/pixel and a FOV of 3.4' x 3.4'. The detector has numerous dead 
and hot pixels. We therefore carefully selected two "clean" areas on the detector for both 
the target and the reference star. TIFKAM also has a known residual charge problem that 
each new image typically has a residual of 0.5-2% of the previous signal. The residual can 
be reduced to <<1% of the original signal when reading out the array several times with 
minimum exposure (see the TIFKAM manual). Therefore, for each science exposure in our 
observations, we readout the detector 3 more times with minimum exposure of 4.29sec and 
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discard the 3 residual frames. The resulting duty cycles are roughly about 68%, 58% and 
46% for 90, 60, and 30 seconds exposures respectively. 

Our first observation started on UT 2010 November 26 at 05:49:55, about 19 minutes 
after the mid-eclipse, and ended at UT 09:10:41, about 134 minutes after the egress. A 
loss of telescope pointing occurred at UT 06:06:00, causing a gap of 17.63 minutes to the 
observation. We started with 60sec exposures and then reduced to 40sec to keep the flux 
within the linearity regime as the airmass decreased. We kept the telescope in focus and 
stayed on the target throughout the observation to minimize instrument systematics. A 
total of 118 images were recorded. The guiding precision was better than 3 pixels for both 
R.A. and Declination. 

The second observation started on UT 2010 November 27 at 06:17:50, about 2 minutes 
after the ingress, and ended on UT 11:38:04, about 149 minutes after the egress. The 
observations on both nights started after the beginning of ingress due to the late rise of 
WASP-12 and technical problems such as pointing & guiding calibration, target centering, 
etc. We started with 90sec exposures for the second night and then reduced to 30sec as 
airmass decreased. A total of 235 images were obtained. The telescope was also kept in 
focus. The guiding drift was about 3 pixels in Declination and 6 pixels in R.A.. The average 
seeing was ~ 1.5" for both nights. 

Reduction of the WASP-12 data followed the same procedures described above for 
CoRoT-1. Due to the sparse field of WASP-12, we only found one good reference star 
(2MASS J06303188+2942273) in the FOV. Nonetheless, the similar spectral type and K 
magnitude of the reference to WASP-12, and the lack of correlation between their fluxes and 
centroid positions on the detector (see §3.2[) still permit reliable flux calibration. Changing 
the exposure times during the observations can also introduce additional systematics to the 
data. However, since the flux of the target star is calibrated with the reference star, both 
stars have exactly the same exposure times, and the flat field and dark frames are well de- 
termined by averaging hundreds of images together, seeing fluctuations are in principle the 
dominating systematics (although other systematics may also exist). Substantial amount of 
"red noise" is thus relatively less likely to be introduced with different exposure times (see 
§3.2p . Therefore, the most significant effect of different exposure times is larger scatter of 
the data points with longer exposures due to more seeing fluctuation^. 

For aperture photometry, we experimented with 48 aperture sizes with a step of 0.5 
pixels. An aperture with radius of 19 pixels for 2010 November 26, and radius of 17 pixels 
for 2010 November 27 presents the smallest scatters in their out-of-eclipse and in-eclipse 



2 This effect can be seen at the beginning of the two observations (see t j3.2[) . 
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Fig. 2. — Flux of CoRoT-1 as a function of X and Y positions of its centroid. The top panel 
shows the correlation of flux with X positions, while the bottom panel shows the correlation 
with Y positions. The linear correlation coefficient is 0.3 for the top panel and for the 
bottom panel. 



data, and thus is used as the final aperture in the photometry. We have also tested the 
eclipse depths with different aperture sizes for both nights in later analysis. Due to the good 
flat fielding of the data and no correlation of flux with centroid positions on the detector, 
eclipse depths are stable and consistent for radii within [-3, +5] pixels of the best aperture 
for both nights. A sky annulus with 29-pixel inner radius and 20-pixel width is used for 
background estimation. Different annulus ranges and sizes have also been explored, and 
indicate consist ent results. The ti mestamp (in UTC) of each image was first converted to 
BJDf/rc- Since ICroll et al.l ( 1201 ll ) has found no evidence of precession for WASP-12b, we 



thus c alculated the orbital phases based on the non-precession ephemerides of ICampo et al. 
(boilh . i.e., period=1.091424days, and transit epoch T ( J BJ£> C7TC )=2454508.97686. 
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Fig. 3. — Final corrected light curve of CoRoT-lb. The first panel shows the X-position 
de-correlated light curve of CoRoT-lb. The second panel shows the residual of the best-fit. 
The bottom panel shows the 15-points averaged light curve. The solid blue lines indicate 
the best-fit light curve of CoRoT-lb. Error bars of the points are calculated from the scatter 
of the data used for averaging. The 3-a percentile of the eclipse depth is between 0.033% - 
0.235%. 



3. Analysis and results 
3.1. CoRoT-lb 

After normalizing the time series of CoRoT-1 with the reference light curve, we still 
see large contaminations of correlated systematics ("red noise") in the data. As can be 
seen in Figure [U the large-scale structures of the time series are highly correlated with the 
centroid positions of the star, possibly caused by a combination of inter-pixel fluctuations 
due to imperfect flat-fielding and other systematic factors, which cannot be corrected by 
the reference light curve since the individual fluctuations of the reference stars are averaged 
out. Figure [2] shows the correlations of CoRoT-l's flux with the X and Y positions of its 
centroid. A clear trend and correlation is visible for the X positions (top panel), although 
the correlation coefficient is low due to the low signal-to-noise of the data. The Y positions 
do not appear to have strong correlation with the flux. 

To correct for these systematic trends, we apply a simple de-correlation function of: 
/ = 1 + a\X, together with a linear slope of the form: 1 + a^i, where / is the flux, X is 



- 9 - 



Before correction 




Fig. 4. — Comparison of CoRoT-lb's noise level with Gaussian expectation. The top panel 
shows the standard deviation of the data as a function of binned points before the X-position 
de-correlation, while the bottom shows the standard deviation after the de-correlation. The 
solid lines indicate the noise levels of the actual data, while the dashed and dot-dashed lines 
indicate the Gaussian noise expectation. 



the positions of the centroid, and t is the orbital phase of each measurement. The term for 
Y positions is not included due to the fact it has little or no correlation with the flux %. 
Because of the insufficient out-of-eclipse baselines, we fit all the in-eclipse and out-of-eclipse 
data together. To account for the eclipse signature, we also fit simultaneously a light curve 
wi th the systemat i c corr ections. The light curve is generated following the prescription 
of iMandel fc Agoll (120021 ) . assuming uniform bodies without limb-darkening. The stellar 
and planetary pa rameters for t he lig ht curve (R p , R s tar, inclination, and semimajor axis) 
are adopted from iGillon et al.l (120091 ) . The free parameters in the least-square fit are: the 
eclipse depth, the level of the out-of-eclipse baseline, and the coefficients ai&a 2 . The known 
durations of ingress and egress are maintained in the fit. The center-of-eclipse timing is also 
fixed to phase=0.50 due to the fact that the orbital eccentricity of CoRoT-lb is consistent 



3 We have also investigated the correction term for Y with the Bayesian Information Criterion and con- 
firmed that the de-correlation function does not prefer such a term for Y positions 
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with zero ( IRogers et al.ll2009t iDeming et al.ll201ll ). and our data have insufficient pre- ingress 
baseline for a robust constraint of ecosa;. 



We employed the Levenberg-Marquardt (LM ) algorithm (IPress et al.lll992l ) for the least- 
square fit. To ensure that we find the global minimum instead of local minima, we searched 
the parameter space extensively with a fine grid of starting points on top of the least-square 
fit. The grid has a few hundred steps for each parameter. The fact that most starting values 
on the grid converge to the same minimum suggests that we indeed have found the global 
minimum. The data points are uniformly weighted such that the xt is nearly 1.0. The global 
best-fit light curve gives an eclipse depth of 0.145%±0.028%. The best-fit model is shown 
by the solid line in the second panel of Figure [H The X-position de-correlated data, the 
residuals of the best-fit and the averaged data are shown in Figure [3J Figure @] compares 
the noise level of CoRoT-lb with the Gaussian noise expectation before and after the X 
de-correlation. The in-eclipse data before the correction (top panel) suggest excessive "red 
noise" above the Gaussian expectation. However, the noise level reduced significantly after 
the correction and becomes more consistent with the out-of-eclipse data that are not affected 
much by the centroid drift (see Figured]). The Gaussian noise level (black dashed line) of 
the in-eclipse data is also improved after the correction. The out-of-eclipse noise level (blue 
solid line), although is only slightly improved, becomes more stable than before at larger bin 
numbers. Overall, the de-correlation has indeed reduced the systematics significantly. The 
average photon noise limit of a single exposure is about 0.049%. Our final precision is thus 
about 4 times of the photon noise limit. 

To verify if the X-position de-co rrelation fun c tion is truly pre ferred, we apply the 
Bayesian Information Criterion^ (BIC)( Liddle 2007 ; Croll et al. 2011 ) for models with and 
without the de-correlation function. The result indicates that the model with the de- 
correlation function has a lower BIC value of 341.3, while the model without de-correlation 
has a higher BIC value of 361.2, suggesting that the X-position de-correlation function is 
indeed superior and preferred. 

We have also investigated the possibility of a quadratic term for the out-of-eclipse base- 
lines using BIC. The results are dependent on the baselines and do not completely justify 
a quadratic term§. In addition, the pre-ingress baseline is also too short to allow a reliable 



4 BIC has been widely used for model identification and selection. Reduces in \ 2 or maximum likelihood 
are penalized for the number of free parameters in BIC. Thus, the model with lower BIC value is generally 
preferred. 

5 We have conducted two BIC tests, one with the post-egress baseline only, and the other with both the pre- 
ingress and post-egress baselines. The first test definitely prefers a linear model for the post-egress baseline 
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quadratic model. Thus, we prefer a linear baseline model instead to avoid possible erroneous 
corrections an inaccurate quadratic term may introduce. A linear baseline can be more reli- 
ably determined since the in-eclipse data have also been used in the joint fit. We also note 
here that a quadratic term, if exists, should only affect the data slightly in large time scale; 
and the time correlated systematics caused by the quadratic term are taken into account in 
the following error analysis. 

To examine the statistical significance and robustness of the eclipse depth and to es- 
timate its error, we conduct 2 statistical tests. We first apply the standard bootstrapping 



technique ( IPress et al.lll992l ). In each bootstrapping iteration, we uniformly resample the 
data with replacement. Typically, a ~ 37% of the original data points are randomly dupli- 
cated in each sample. For each new sample, we re-fit the X de-correlation function and the 
linear slope simultaneously with the light curve to determine the eclipse depth, using the 
aforementioned grid search and LM minimization. This technique is suitable for unknown 
distributions like our case, and can robustly test the de-correlation model and the distribu- 
tion of the parameters. A total number of 1500 iterations are performed and the resulting 
distribution of the eclipse depth is nearly Gaussian, with a median and l-a deviation of 
0.146%±0.027%, highly consistent with the previous best-fit. 



For the second test, we use the "prayer-bead" residual permutation method (I Winn et al. 



20081 . and references therein). In brief, we subtract the best-fit model from the data and 
shift the residuals pixel-by-pixel. The shifted residuals are then added back to the best-fit 
model to simulate a new set of data. The same de-correlation function and light curve are 
then employed to re-fit the new data for each iteration. We also reverse the residuals and 
iterate this process again, resulting in a total number of 751 iterations (i.e., 2N-1, where 
N=376 is the number of data points). This method maintains the time-correlated errors and 
is therefore another robust way of testing our fit. Due to the un-corrected "red noise" in 
the residual, the eclipse depth shows larger scatter in this test, and the final distribution is 
top-flat. The resulting median depth and l-a error is 0.148%±0.049%, and the 3-cr percentile 
is from 0.033% to 0.235%, suggesting the eclipse depth is detected at 3-cr significance and is 
consistent with previous results. The results of these analysis are summarized in Table [TJ 

As a final test and cross check of our eclipse signal, we conduct a least-square fit to 
the original data without the X-position decorrelation but only a light curve and a linear 



(i.e., significantly lower BIC value). The second test marginally prefers a quadratic baseline. However, the 
quadratic fit is strongly leveraged by a few points (6) at the very end of the observation that have large 
scatters due to deteriorated seeing. If these points are excluded, the BIC once again strongly prefers a linear 
baseline. Therefore, a quadratic baseline cannot be completely justified. (When including these points, the 
quadratic baseline results in a best-fit eclipse depth of 0.102±0.033%, still consistent with our final result.) 
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slope. The resultant best-fit depth is 0.125%±0.029%, while the result from bootstrapping 
is 0.127%±0.028%, both are consistent with the previous results within error bars. The 
residual permutation method is also applied in this test, leading to a depth of 0.121% with 
a 3-cj percentile between 0.019% - 0.029% - still consistent with the other tests. The larger 
error in this permutation test is expected because of the excessive un-corrected "red noise" . 

Based on the above tests and their consistency, we conclude that our detection of the 
eclipse of CoRoT-lb is real and has at least 3-cr significance. We report the final eclipse 
depth as 0.145%±0.049% based on the result from the original best-fit and the largest error 
from the residual permutation test (see Tabled]). 



3.2. WASP-12b 

The two light curves of WASP-12b are analyzed in a similar way as CoRoT-lb. The 
reduced light curves of WASP-12b after normalizing with the reference star are shown in 
the top panels of Figure [51 To correct for the systematics, we first investigate if centroid 
de-correlation functions are preferred by BIC. The calculations show that models without 
centroid de- correlation have lower BIC values than those with X or Y de- correlations for 
both nights, suggesting that simpler models are preferred. Figure M shows the correlations 
of the reduced flux of WASP-12b (i.e., shown in top panels of Figure E]) with the X and Y 
positions of its centroid. The figures suggest there is no obvious correlations between its 
flux and centroid positions. Therefore, no X or Y de-correlation terms are required in the 
models, as already suggested by the BIC tests. 

To investigate if quadratic terms are preferred in the out-of-eclipse baseline models of 
WASP-12b, we also calculate the BIC values for each night. Since there is no pre-ingress 
baseline for either night, only the post-egress baseline is used. It turns out that for both 
nights, linear background models are strongly preferred to quadratic models based on their 
lower BIC values. Therefore, in our final analysis of the WASP-12b data, we only apply a 
simple linear baseline model together with a light curve. 



Because WASP-12b's orbital eccentricity is consistent with zero (jCroll et al.ll201ll ) and 
we have no pre-ingress baseline, we also keep the center-of-eclipse timing fixed to phase=0.5, 
and maintain the known durations of ingress and e gress in the fit. Stellar parameters for the 



model light curves are adopted from iHebb et al.l (120091 ) . An extensive grid of the starting 



parameters is also applied on top of the LM minimization to ensure we find the global 
minimum. The data points are weighted uniformly such that xt is about 1. The final best- 
fit eclipse depth is 0.281%±0.085% for 2010 November 26, and 0.316%±0.079% for 2010 
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November 27. 

The best-fit models, including the light curves, are shown in the top panels of Figure 
(blue lines). In the two top panels, both nights show a linear systematic trend in the 
background, which can be caused by large time scale seeing variation and/or drift of the 
thermal background. The middle panels of Figure [5] show the background-corrected data 
along with the best-fit light curves. Corresponding residuals and the averaged data are 
shown in the third and the bottom panels. Figure [7| compares the noise level of WASP-12b 
with the Gaussian noise expectation before and after the background correction. Overall, 
there is very low systematic "red noise" in both nights, and their standard deviations follow 
the Gaussian expectations closely. The bottom panel of 2010 November 26 shows that after 
the background correction, the noise level is effectively reduced to the Gaussian expectation. 
For the 2010 November 27 data, although the noise level stays nearly the same, the binned 
standard deviation becomes more stable at large bin numbers after the correction. 

To examine the statistical significance and robustness of the model, we also implement 
the same bootstrapping and residual permutation tests for WASP-12b. A total number 
of 1500 iterations are carried out for the bootstrapping test, resulting in a Gaussian-like 
distribution with a median depth and 1-a deviation of 0.290%±0.085% for 2010 November 
26 and 0.329%±0.077% for 2010 November 27. The residual permutation tests have 235 
and 469 iterations for November 26 and November 27 respectively, resulting in a depth of 
0.279%±0.077% for the former and 0.318%±0.095% for the latter. The results are summa- 
rized in Table [TJ The eclipse depths from both nights and all three methods are consistent 
with each other and all suggest better than 3-cr significance. In addition, since the system- 
atics of the two nights (e.g., airmass, seeing variation, thermal background, pointing drifts, 
observation timing and duration, etc.) are very different, and the only common feature of 
them is the eclipse signal, we conclude that the above detections are real. Finally, we use the 
eclipse depths from the least-square fits and the largest error of the 3 methods to combine the 
two nights together. The final joint e clipse depth of WA SP-12b is 0.299%±0.064% (4.7-a), 



consistent with the Ks band result of iCroll et all f l201lh . 0.309%+J^J. 



4. Discussion 



We summarize our new H band measurement of CoRoT-lb in Figure [H] along with 
other previous results. Our H band measurement corresponds to a brightness temperature 
of 2280^23oK for CoRoT-lb, slightly higher than its equilibrium temperature of 2180K, as- 
suming a zero Bond Albedo and no heat redist ribution from the da yside to the nightside. The 
temperature-inverted atmospheric models of iGillon et al.l (120091 ) and iDeming et al.l (120111 ) 
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are fu lly consistent with our result. Our measurement supports the conclusion of lGillon et al. 



( 120091 ) that a model with uniform redistribution of stellar flux across the entire planet surface 
is too cool to match the data (the dashed line). We also find a blackbody of T=2380±100K 
with no heat redistribution best fits the data at all wavelengths, implying an isothermal 
region across much of the photosphere and very ine fficient transport of heat to the night 
side. This is also in agreement with the conclusions of iDeming et all ( 1201 ll ). The addition of 
the new H band measurement is not sufficient to differentiate the best-fit models of previous 
studies at other wavelengths. Thus, in order to put more stringent constraints to current 
models, measurements with better precision at these or new wavelengths are necessary. 



5. Conclusions 

We have made a new detection of the very hot Jupiter CoRoT-lb's thermal emission at 
the H band. The 3-a detection suggests a eclipse depth of 0.145%±0.049%. This result is 
consistent with the conclusions of previous studies that the planet probably has a thermal 
inversion layer at high altitude, and has an isothermal region with inefficient heat transport 
across its dayside and nightside. 

We have also detected the thermal emission of WASP-12b at the Ks band on two nights, 
at a 4.7-cr joint ecli pse depth of 0.299 %±0.065%. This result independently confirms the 



previous detection of lCroll et all (120111 ) with a different telescope and instrument, suggesting 
the robustness of both measurements and also validating our data reduction and analysis 
method. 

Although more precise measurements are still required to better constrain models for 
CoRoT-lb and WASP-12b, our detections of the two planets' secondary eclipses have brought 
the Palomar 200-in telescope and the MDM 2.4m telescope to the inventory of telescopes 
that have demonstrated the capability of detecting hot Jupiter's thermal emission from the 
ground. Together with previous results made with other telescopes, these detections suggest 
that ground-based observations are now mature and becoming a widely available tool to 
characterize the atmospheres of hot Jupiters. 
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Fig. 5. — Light curves of WASP-12b. The top panels show the flux of WASP-12b normalized 
to the reference star, overplotted with the best-fit background correction models and light 
curves. The middle panels show the background-corrected flux with the best-fit light curves. 
Residuals of the best-fit are shown in the third panels. The averaged data are shown in the 
bottom. Error bars are calculated from the scatter of the data used for averaging. 
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Fig. 6. — Flux of WASP-12b as a function of X and Y positions of its centroid. The left two 
panels show the correlation for 2010 November 26, while the right panels show the correlation 
for 2010 November 27. Y positions (top panels) correspond to R.A. for MDM/TIFKAM, 
while X positions (bottom panels) correspond to Declination. 



Table 1. CoRoT-lb and WASP-12b eclipse depth and error estimate 



Target 


Method 


Eclipse Depth 3-cr percentile 




Least-square fit 


0.145%±0.028% 0.061% - 0.229% 


CoRoT -lb 


Bootstrap 


0.146%±0.027% 0.065% - 0.227% 




Residual Permutation 


0.148%±0.049% 0.033% - 0.235% 




Final result 


0.145%± 0.049% 




Least-square fit 


0.281%±0.085% 0.025% - 0.535% 


WASP-12b 


Bootstrap 


0.290%±0.085% 0.035% - 0.545% 


(2010 November 26) 


Residual Permutation 


0.279%±0.077% 0.124% - 0.422% 




Final result 


0.281%± 0.085% 




Least-square fit 


0.316%±0.079% 0.079% - 0.533% 


WASP-12b 


Bootstrap 


0.329%±0.077% 0.098% - 0.560% 


(2010 November 27) 


Residual Permutation 


0.318%±0.095% 0.113% - 0.514% 




Final result 


0.316%± 0.095% 


WASP-12b 


Joint solution 


0.299%± 0.064% 
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Fig. 7. — Comparison of WASP-12b's noise level with Gaussian noise expectation. The top 
figure shows the standard deviations of the out-of-eclipse data as a function of binned points 
for 2010 November 26. The bottom figure shows the standard deviations of both the in- 
eclipse data and the out-of-eclipse data for 2010 November 27. The solid lines indicate the 
binned down noise levels of the actual data, while the dot-dashed lines indicate the Gaussian 
noise expectation. The in-eclipse data of 2010 November 26 are not shown due to the small 
number of points. 
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Fig. 8. — Comparison of atmospheric models with CoRoT-lb data. The new H band data is 
shown in orange. Previous measurements in other bands are shown in blue. The black line 



i ndica tes the best-fit blackbody model. The temperature-invert ed model fromlDeming et al. 
( 1201 ll ) is shown by the green line. The dayside only model from lGillon et al. is shown 

by the solid purple line, while the model with full heat distribution over the day and night 
sides is shown in dashed purple line. The inverted transmission profiles of the measurements 
are shown on the top of the figure. The dotted line that intersects the CoRoT b a ndpas s on 
the top left corner indicates the cut-off of the blue channel used in lSnellen et al.l ( 2009 ). 



